Critical behavior of Binder ratios and ratios of higher order cumulants 
of conserved charges in QCD deconfinement phase transition 
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Binder liked ratios of baryon number are firstly suggested in relativistic heavy ion collisions. Using 
3D-Ising model, the critical behavior of Binder ratios and ratios of higher order cumulants of order 
parameter are fully presented. Binder ratio is shown to be a step function of temperature. The 
critical point is the intersection of the ratios of different system sizes between two platforms. From 
low to high temperature through the critical point, the ratios of third order cumulants change their 
values from negative to positive in a valley shape, and ratios of fourth order cumulants oscillate 
around zero. The normalized ratios, like the Skewness and Kurtosis, do not diverge with correlation 
length, in contrary with corresponding cumulants. Applications of these characters in search critical 
point in relativistic heavy ion collisions are discussed. 
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One of the main goals of current relativistic heavy ion 
experiments is to locate the critical point of QCD de- 
confinement phase transition. The critical character is 
that the correlation length £ goes to infinite larger at 
infinite system. For finite system, like the formed one 
in relativistic heavy ion collisions, the correlation length 
should be a finite maximum. Therefore, the various cor- 
relation length related observables are suggested in rela- 
tivistic heavy ion collisions [l|. 

It has been recently shown that near the critical point, 
the density-density correlator of baryon-number follow 
the same power law behavior as the correlator of the 
sigma field, which is associated with the chiral order pa- 
rameter 0, Q . Therefore, the baryon number is consid- 
ered as an equivalent order parameter of formed system 
in nuclear collisions |4|. 

From statistic physics, it also shows that the suscepti- 
bilities of order parameter is directly related to the fluc- 
tuations of conserved charges, i.e., 



(SN*) = VT Xi . 



(1) 



Xi is the ith order susceptibility. (SN 1 ) = ((N — N) 1 ) 
is the ith order cumulants of the conserved charge num- 
ber N. For three flavor QCD, the conserved charges are 
baryon-number, strangeness, and electric charge [5J. 

The third and forth order cumulants of conserved 
charges are defined respectively as, 



K 3 = (SN 3 ), K 4 = (8N A ) - 3(6N 2 ) 2 



(2) 



In the vicinity of critical point, they are argued to be pro- 
portional to the higher power of correlation length, i.e., 
£ 4 - 5 and £ 7 [HQ, respectively. So they are more sensitive 
to the correlation length, and highly recommended. 

In experiments [8j], properly normalized cumulants, 
i.e., Skewness and Kurtosis, 
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are actually presented. As the second cumulant is also 
proportional to a certain power of correlation length Q , 
if such normalized Skewness and Kurtosis diverge with 
correlation length is not clear from the theoretical point 
of view. 

From theoretical side, the ratios of high order cumu- 
lants to the second one, e.g., 



R 3 ,: 



(SN 3 ) 
(SN 2 )' 
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3(5N 2 ). 
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are estimated |10l - ll4| . The Lattice QCD with two light 
quark degrees of freedom shows that these ratios of the 
baryon number, strangeness, and electric charge have 
pronounced peaks from low to high temperature in the 
transition region of chiral symmetry break [loj . The 
effective models in the mean-field approximation also 
shows that there are peak, valley, and oscillating struc- 
tures near the deconfinement and chiral phase transi- 
tions [U EH [Hj]. However, all these are obtained un- 
der some approximations due to the difficulties in Lattice 
QCD calculations 15] and model estimations (l3j . 

Although the concrete form of interactions varies from 
one system to another, according to the theory of univer- 
sity, the critical exponent of equivalent measurement is 
identical in the same university class. This allows us to 
study the critical behavior of complex system by known 
simple one. 

It is known that the QCD deconfinement phase tran- 
sition corresponds to the restoration of 0(1) symmetry, 
which is the same university class of 3D-Ising model [1 6j j . 
Therefore, the critical behavior of all above mentioned 
cumulants of baryon number can be easily obtained from 
the corresponding cumulants of order parameter in 3D- 
Ising model. 

Moreover, it it known in statistical physics that the 
Binder ratio of order parameter is a direct location of 
critical point [l7j . Generally, the Binder liked ratios are 
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normalized raw moments of order parameter. The third 
and fourth Binder liked cumulant ratios can be simply 
defined as, 

Here we take 3D-Ising model as an example. The order 
parameter in the model is the magnet M = Y2i=i sI/Nl 
of spin s in all lattice sites Nl . 

Equivalently, the order parameter in relativistic heavy 
ion collisions is the baryon number. The temperature, 
or the controlling parameter, is the incident energy. The 
size of the formed system is mainly determined by the 
overlapped area, i.e., centralities. So if we pass through 
the region of critical incident energies in relativistic heavy 
ion collisions, the Binder ratios of baryon number can 
be served as a good location of critical point of QCD 
deconfmement phase transition. 

In this paper, we firstly present the critical behavior of 
Binder ratios in 3D-Ising model, and demonstrate why 
they are helpful, in particular, in locating the critical 
point in relativistic heavy ion collisions. Then, the crit- 
ical behavior of Skewness, Kurtosis, i?3,2, and R^2 are 
presented and discussed, respectively. Meanwhile, from 
finite-size scaling of the susceptibilities, the critical be- 
havior of those ratios are estimated model independently. 
Finally, the conclusions are drawn. 

The critical behavior of Binder ratios, B3 and B4, 
in 3D-Ising for 4 different lattice sizes are presented in 
Fig. 1(a) and (b), respectively. Where the simulation of 
3D-Ising model is based on the wolff algorithm We 
can see that both B3 and B4 show a step jump in the 
vicinity of critical temperature. The physical meaning of 
this jump is clear. 

When the temperature is much lower than the critical 
one, the system is almost order and the fluctuation of 
order parameter is very small, i.e., 

(M n ) ~ (M) n (for n = 2,3,4- ••)■ (6) 

So it results the lower platform, which is 1 for all orders 
of Binder ratios at all system sizes, as shown in Fig. 1. 
When the temperature approaches to critical one, the 
correlation length starts to increase with temperature 
and the fluctuations become larger and larger. Their 
critical behavior is system size dependent and described 
by finite-size scaling. Only at critical temperature, all 
size curves intersect to the fixed point, where they are 
system size independent [19], as shown in Fig. 1. When 
the temperature is much higher than the critical one, the 
system is totally disordered. It approaches again to a 
constant. This forms the platform at high temperature. 
It is 1.6 and 3 times larger than the lower platforms for 
the third and fourth order Binder ratios, respectively. So 
the higher the order of Binder ratio, the larger the gap 
of the step function. 

This step function liked behavior can be served as a 
very good probe of critical point in relativistic heavy ion 
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FIG. 1: (Color online) The temperature dependence of Binder 
ratios in Eq. ([5| in the vicinity of critical temperature in 3D- 
Ising model for 4 different lattice sizes. 
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FIG. 2: (Color online) The temperature dependence of Kur- 
tosis (a) and Skewness (b) in Eq. (|3} in the vicinity of critical 
temperature in 3D Ising model for 4 different lattice sizes. 

collisions, where critical incident energy is difficult to as- 
sign precisely in priori. So if we scan incident energies, 
and observe two platforms at low and high energy re- 
gions, respectively, then the critical one is most probably 
between them. We can finely tune the incident energy 
in the region and precisely determine the critical energy 
and exponents. 

The Skewness and Kurtosis of order parameter in 3D- 
Ising model for 4 different lattice sizes are presented in 
Fig. 2(a) and (b), respectively. We can see from the fig- 
ure that they change sharply in the vicinity of the critical 
temperature. The Skewness first drops down and then 
goes up, and Kurtosis oscillates with temperature. Their 
values are system size dependent. Their signs change re- 
spectively near the critical point. Former in Fig. 2(a) 
changes from negative to positive when the temperature 
is increased through the critical point, while the later in 
Fig. 2(b) becomes negative only when the temperature 
is close to the critical point. The sign change in Skew- 
ness, or third order cumulants, is expected in effective 
models pHHHl- 

As we know that the Skewness and Kurtosis measure 
the symmetry and sharpness of the distribution, respec- 
tively. The distributions of order parameter M near 
the critical point at system size L = 8 are shown in 



3 



0.06 




FIG. 3: (Color online) The distributions of order parameter 
near critical temperatures in 3D-Ising model at system size 
L = 8. 
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FIG. 4: (Color online) The temperature dependence of i?3,2 
(a), and ^4,2 (b) in the vicinity of critical temperature in 
3D-Ising model for 4 different lattice sizes. 



Eq. © are, 



Fig. 3. Where we can clearly see that the long tail of 
the distributions changes from the left to the right side 
when the temperature is increased through the critical 
point, and the peak of the distribution vary from sharp 
to flat when temperature is approached the critical point. 
The same trend has been observed in percolation model, 
in studying clusterization phenomena in nuclear multi- 
fragmentation (2lj . 

This character can also be served as a signal associ- 
ated with the appearance of critical point in relativistic 
heavy ion collisions. If we observe sign change of Skew- 
ness (Kurtosis) of baryon number at a certain incident 
energy region, it most probably predicts the appearance 
of critical point in the nearby incident energy region. 

The Skewness and Kurtosis also converge to two con- 
stants when the temperature is away from critical point, 
as shown in Fig. 2(a) and (b). But the constants at low 
and high temperatures are close to zero and 1, respec- 
tively. The gap between them are small and does not 
change very much with the order of cumulants, unlike 
the Binder ratio. 

Moreover, all size curves of Skewness (Kurtosis) inter- 
sect at critical point. This can be easily understood from 
finite-size scaling of susceptibilities, i.e., 



Xi (t,L)=L^ v P Xi (tL^)- 



(7) 



Where the 7$ is the critical exponents of ith order sus- 
ceptibility, and v is the critical exponent of correlation 



length = t v at infinite system, t 



T-T c 



is reduced 



temperature, and T c is critical temperature. In the vicin- 
ity of critical point, the correlation length at finite system 
is approximately the same order of the system size, i.e., 

For X3 an d Xii 1t>I v = 4-5, 74/^ = 7, respectively @. 
So the critical behavior of the Skewness and Kurtosis in 



Ki/Kl 



VT X3 



L 3 L^P X3 (tL^) 



{VTf/ 2 xl /2 L*- 5 L 3 TV 2 P x l 2 {tL 1 /") 
T-^FsitL 1 ^), 
VT Xi 



L 3 L 7 P x JtL 1 / 1 ') 



{VTfxl 
T^FKitL 1 /"). 



L*L 4 TPl 2 



(8) 



They no long diverge with correlation length, or system 
size. At the critical temperature t = 0, the scaling func- 
tion, i.e., Fs(0) or Fk(0), is system size independent con- 
stant. All size curves intersect to the constant, i.e., the 
fixed point [19j . 

From this simple estimation and Fig. 2, we can see 
that normalized high order cumulants, i.e., Skewness and 
Kurtosis, do not directly diverge with correlation length 
any more, different from corresponding cumulants, K3 
and K4, which are proportional to £ 4 ' 5 and £ 7 , respec- 
tivelyJsQ. 

The i?3,2, and i?4,2 of order parameter in 3D-Ising 
model for 4 different lattice sizes are presented in 
Fig. 4(a) and (b), respectively. We can see again from 
Fig. 4(a) that R^^ changes its value sharply from neg- 
ative to positive when temperature is increased through 
the critical point. R^^ in Fig. 4(b) oscillates greatly with 
temperature near the critical point. These qualitative 
features, i.e., sign change in third moment, and oscil- 
lating structure in forth cumulants, are consistent with 
estimations of effective models [HI, Q, H(| . 

i?3_2 and R4.2 are very sensitive to the system size, or 
correlation length. Their values become very large when 
system size increases. The critical exponent of -R3.2, and 
i?4 ; 2 can be roughly estimated from finite-size scaling of 
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susceptibilities, i.e., 

VT X3 = L^P X3 {tL 1/v ) 
3,2 VT X2 L*?P Xt (tLV») 

VT X i I?eP Xi {tL l ' v ) 
4 ' 2 VT X2 L*?Pn(tlM») 

= eF^tiM"). (9) 

So R^ : 2 and Ra,2 diverge with correlation length as £ 2 5 
and £ 5 , respectively. 

In this paper, the measurements of Binder liked ratios 
of conserved charges are firstly suggested in relativistic 
heavy ion collisions. Using 3D-Ising model, it is shown 
that near the critical temperature, Binder ratios is a step 
function of temperature. The gap of the step function 
is 1.6 and 3 times wider for the third and forth order 
Binder ratios, respectively. This can be served as a good 
identification of critical behavior in relativistic heavy ion 
collisions, where the critical incident energy is unknown 
in prior. The critical point is the intersection of Binder 
ratios of different size systems between two platforms. 

The critical behavior of Skewness, Kurtosis, R^^ and 



i?4 ; 2 at various system sizes are also studied by 3D-Ising 
model, and estimated by finite size scaling. When the 
temperature is increased through the critical point, the 
ratios of the third order cumulants change their values 
from negative to positive in a valley shape, and ratios of 
fourth order cumulants oscillate around zero. All size 
curves of Skewness (Kurtosis) intersect at the critical 
point. The normalized ratios, like the Skewness and Kur- 
tosis, do not diverge with correlation length. While, un- 
normalized ratios, ^3,2 and -R4.2, are divergent with cor- 
relation length. They are proportional to £ 25 and £ 5 , 
respectively, and very sensitive to the system size near 
the critical temperature. 

These critical characters may show up at the energy 
dependence of corresponding ratios of conserved charges. 
Their behavior at coming relativistic heavy ion experi- 
ments at RHIC, SPS, and FAIR are called for. 
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